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Abstract 

The one-dimensional deterministic economic model recently studied by Gonzalez- 
Estevez et al. [Physica A 387, 4367 (2008)] is considered on a two-dimensional 
square lattice with periodic boundary conditions. In this model, the evolution of 
each agent is described by a map coupled with its nearest neighbors. The map has 
two factors: a linear term that accounts for the agent's own tendency to grow and 
an exponential term that saturates this growth through the control effect of the 
environment. The regions in the parameter space where the system displays Pareto 
and Boltzmann-Gibbs statistics are calculated for the cases of von Neumann and of 
Moore's neighborhoods. It is found that, even when the parameters in the system 
are kept fixed, a transition from Pareto to Boltzmann-Gibbs behavior can occur 
when the number of neighbors of each agent increases. 
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1 Introduction 



In the last few years, different probabilistic models [T][2][5|HI5] have been pro- 
posed to explain the wealth distribution in western societies pmiSMTUfTT] . 
namely the Boltzmann-Gibbs distribution grouping about the 95% of individ- 
uals, corresponding to those belonging to the low and middle economic classes, 
and the Pareto distribution consisting of the 5% of individuals possessing the 
highest wealths. The majority of these models explain the wealth distribution 
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as the consequence of random processes (TO]. However, some misconceptions 
lie along this line of research, since most people are aware that almost every 
economic trade is done under the rationale force of some interest or some final 
profit, and only occasionally this is conducted by chance. Thus, in order to 
shed light on the problem of how the wealth distribut! es in human society, 
it is important to dispose of other economic models incorporating different 
degrees of determinism in the interaction among the agents. 

One of these models was recently proposed [12] and studied in detail for the 
one-dimensional case in [13]. This is a completely deterministic model that 
reproduces realistic wealth distributions, i.e., the Pareto and the Boltzmann- 
Gibbs (BG) distributions, for different values of parameters. Moreover, it is 
possible to produce a transition from Pareto to BG behavior by only modi- 
fying one parameter. It means that it is possible to bring the system from a 
particular wealth distribution to another one with a lower inequality by per- 
forming only a small change in the system configuration. This is an advantage 
respect to other random models where it is necessary to perform a structural 
reconfiguration of the system in order to get this type of transition [2][3] . 

In this paper, we ask if some other strategies can be implemented in this model 
[T^fTB] in order to induce a transition in its asymptotic statistical behavior. 
We proceed to affirmatively answer this question by showing that different lo- 
cal groupings of the agents that occur, for instance, by changing the topology 
of the lattice or the number of neighbors of each agent, can lead to a transi- 
tion from the Pareto to BG behavior even when the parameters of the system 
are kept fixed. Hence, the first step we undertake in Section 2 is the charac- 
terization of the statistical behaviors that this model exhibits on its space of 
parameters when it is implemented on a two-dimensional lattice. Particularly, 
the cases of von Neumann and Moore neighborhoods are studied in detail. 
Then, in Section 3, the results obtained for these configurations are compared 
with those previously found for the one- dimensional case. The regions on the 
parameter space where a transition between different statistical behaviors can 
take place, are identified in this section. Section 4 contains our conclusions. 



2 The deterministic economic model on a two-dimensional lattice 

The system [T2][T3] consists of agents placed at the nodes of a network. Here, 
the network is a square lattice with periodic boundary conditions. Each agent, 
representing an individual, a company, a country or other economic entity, is 
identified by a pair of indexes with i,j = 1, . . . , N. The dynamics of 

each agent is described by a discrete-time map that expresses the competition 
between its own tendency to grow and an environmental influence that controls 
this growth. The dynamics of the system is described by the coupled map 
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where x]'-' > gives the state of the agent at discrete time t, and it may 
denote the wealth of this agent; the factor Vijxl'-' expresses the self- growth 
capacity of agent characterized by a parameter rjj-; represents the 

local field acting at the site at time t; lyiijj) is the set of agents in the 
network coupled to agent and r]{i,j) is the cardinality of this set; and 

ttij measures the coupling of agent (z, j) with its neighborhood; it can also be 
interpreted as the local environmental pressure exerted on agent {i,j) . The 
negative exponential function acts as a control factor that limits this growth 
with respect to the local field. With the dynamics given by Eqs. ([H) the largest 
possibility of growth for agent {i,j) is obtained when xl'-' ~ ajj\E'J'"', i.e., when 
the agent has reached some kind of adaptatio! n to its local environment. 

We consider two different nearest-neighbor interactions: the 4-cell von Neu- 
mann neighborhood, 

= y^t + + xi' + xi' j ; (2) 



and the 8-cell Moore neighborhood, 
1 

t ^\ t t t t 

j-l,j-l I j-l,j + l I j+l,j-l I j+l,j + lN 



We focus on a homogeneous system where all agents possess the same growth 
capacity, rjj- = r, and are subject to a uniform selection pressure from their 
environment, = a. Thus, the parameter a can be interpreted as a homo- 
geneous social constraint on the agents to reach a wealth level proportional 
to that of their environment. The value a = 1 implies a tendency towards 
being totally balanced with the neighborhood. The cases a < 1 corresponds 
to a lower level of expectation among agents that restricts the possibilities 
for improving their relative wealth. When a > 1 the agents possess a greater 
stimulus for overcoming their local neighbors. 

We study the collective behavior of the system described by Eqs.([T]) in the 
space of parameters (a, r). For all the simulations shown the system size is 
316 X 316 (A^ ~ 10^) and the values of the initial conditions are uniformly 
distributed at random in the interval Xq-' G [1, 100]. Also, a transient of 10^ 
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Fig. 1. Left panels: Regions where Boltzmann-Gibbs behavior P{x) ~ ap- 
pears on the space of parameters (a, r), indicated by the label BG, for the (a) 
4-cell and (c) 8-cell neighborhoods cases. The grey code on the right indicates the 
values of the scaling exponent fi. Right panels: Regions where Pareto behavior 
P{x) ~ x~" occurs on the plane {a,r), labeled by P, for the (b) 4-cell and (d) 
8-cell neighborhoods cases. The grey code on the right indicates the values of the 
exponent a. 



iterations is discarded before arriving to the asymptotic regime where all the 
calculations are carried out. When indicated, time averages are done over the 
next 100 iterations after the transient, and this result is newly averaged over 
100 different realizations of the initial conditions with the same parameter 
values. 



Figure 1 shows the regions where the probability distribution P{x) displays 
BG behavior, P{x) ~ e~^^, and Pareto behavior, P{x) ~ x""", in the space of 
parameters (a, r) for the system Eq. ([1]). The parameters a and r are varied in 
intervals of size Aa = 0.02 and Ar = 1, respectively. For each pair (a,r) and 
after undergoing the averaging process described above, semilog and log-log 
linear regressions are used to calculate the exponents /i and a for the corre- 
sponding distributions, and only those results yielding a correlation coefficient 
greater than 0.96 are shown in Fig. 1. 

In comparison with the one-dimensional case [13j, the BG and Pareto regions 
in the space of parameters are larger for the square lattice. In both cases, 
the BG behavior occurs for the lower values of the parameter a, but in two- 
dimensions there is also a narrow BG region for greater values of a. In general, 
when the value of a increases, the population enters in a competitive regime 
that provokes the appearance of the Pareto behavior. The scaling exponents 
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Fig. 2. Left panels: Mean field Ht of the system at t = 10^ for the Boltzmann-Gibbs 
(BG) and Pareto (P) regions for the (a) 4-cell and (c) 8-cell neighborhoods cases. 
The grey code on the right indicates the values taken by Ht. Right panels: The 
quantity (a) on the space of parameters (a, r), for the (b) 4-cell and (d) 8-cell 
neighborhoods cases, with the grey code indicated on the right. 

obtained for both regions are similar to the one-dimensional case. Those ex- 
ponents that correspond to the Pareto behavior are in the range a G [2.4, 3.0], 
in agreement with those found in actual economic data [7|15|16j . 

The mean field of the system or average wealth per agent at a time t is defined 

as 

1 ^ . . 

Ht = j^T.^t'- (4) 



Similarly to the results reported in [13\ for the one-dimensional case. Fig- 
ure 2(a-c) shows the asymptotic value of the global mean field for the regions 
where BG and Pareto behaviors are observed on the space of parameters (a, r), 
for the 4-cell and 8-cell neighborhood cases. Note that, although the values 
of the initial states of the agents are randomly distributed on the interval 
[1,100], the system evolves to an asymptotic state where Ht takes values on 
the smaller interval [0.6,3.25]. Let us observe that the BG region presents a 
higher mean wealth than the Pareto one, telling us that in this model the 
BG behavior would be desirable for the general benefit of the ensemble. In 
the same line, the BG and Pareto regions with 8-cell neighborhood display a 
slightly bigger mean wealth than the 4-cell neighborhood and the ID cases. 
It seems to indicate that the increasing number of the local neighbors of each 
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agent can generate more richness in the system.! On the other hand, for some 
values of the parameters a and r, the states of agents in the system at a given 
time exhibit a large dispersion. For those parameters, the values of the state 
of any agent present large fluctuations over long times. To characterize these 
fluctuations, we define the instantaneous standard deviation of the mean field 
as 



o-t 




1/2 



(5) 



After discarding 10^ transients, we calculate the mean value of at over 100 
iterations, and then average this result over 100 realizations of initial condi- 
tions. The resulting average dispersion, denoted by (a), is shown in Fig. 2(b-d) 
on the plane (a, r). Note that for some regions of parameters the quantity (at) 
can be much greater than the value of the mean field. 

The large dispersions observed in Fig. 2 (b-d) reflect the inequality in the 
wealth distribution among agents in this system. To characterize the degree 
of inequality in the wealth distribution we use the Gini coefficient defined at 
a time t as [T7] 
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Fig. 3. (a) Gini coefficient at t = 10^ as a function of the parameters a and r for the 
(a) 4-cell and (b) 8-cell neighborhoods cases. The labels BG and P indicate the 
regions with Boltzmann-Gibbs and Pareto behaviors. The grey code on the right 
indicates the values taken by the Gini coefficient. 

A perfectly equitable distribution of wealth at time t, where x]''^ = xl'^ V j, fc, /, 
yields a value Gt = 0. The opposite situation, where one agent has the total 

N . . 

wealth J2 xl'-' , has a value of Gt = 1. Figure 3 shows the asymptotic value 

of the Gini coefficient on the plane of parameters {a,r). Note that the Gini 
coefficient reaches larger values, i.e. Gt G [0.65, 0.85] and Gt G [0.7, 1] for 
the 4-cell and 8-cell neighborhoods, respectively, in the regions associated to 
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Pareto regimes, while it takes lower values, i.e. Gt G [0.35,0.65], in the re- 
gion corresponding to BG behavior. This results agree with our qualitative 
understanding that equity is more favored in the presence of a larger middle 
economic class in a society, as expressed by a BG distribution. 



3 Transition from Pareto to Boltzmann-Gibbs behavior 

As it was put in evidence in the one- dimensional realization of this economic 
system studied in [13], it is possible to produce a transition from a Pareto 
to a BG behavior (or viceversa) maintaining its structural configuration and 
performing only a variation in the parameters (a, r) of the system. In the two- 
dimensional it can be seen in Fig. 1, it is also possible to make such 
transition, for example, by changing r with a fixed value of a. We call the 
transition from Pareto to BG behavior (or viceversa) induced by varying the 
parameters the Transition Strategy I. 

Figure 4 reveals another mechanism capable of provoking a change in the 
statistical properties of this economic model. We call it the Transition Strategy 
II. It consists in performing a reconfiguration of the local connectivity of the 
agents. Figure 4(a) shows the region of the parameter space (a, r) where a 
change of the topology of system from a one- dimensional array to a two- 
dimensional square lattice with 4-cell neighborhood takes the system from a 
Pareto to a BG asymptotic distribution. The same transition can be reached 
in the region plotted in Fig. 4(b) by transforming the one-dimensional array 
to a two-dimensional lattice with 8-cell neighborhood. Figure 4(c) displays the 
region of parameter space where a modification of the number of neighbors, 
from a 4-cell von Neumann neighborhood to a 8-cell Moore neighborhood, 
generates a Pareto-BG transition in the two-dimensional array. 

As shown in Fig. 4(d), the transitions in statistical behaviors achieved through 
either type of strategy produce a decrease of the Gini coefficient and, therefore, 
a more equitable wealth distribution in the system. 



4 Conclusions 

Nowadays, both the exponential and power law statistical behaviors in the 
economies of western countries are well documented. Different models for the 
interaction of the economic agents have been reported in the literature. Those 
models with a probabilistic inspiration must apply a strong reconfiguration in 
the structural properties of the system in order to obtain a transition from a 
BG to a Pareto distribution in its asymptotic state. 
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Fig. 4. Regions in the parameter space (a, r) where the deterministic model displays 
a transition from Pareto to Boltzmann-Gibbs statistical behavior through a recon- 
figuration of the local neighborhoods, (a) The system changes its topology from a 
one-dimensional array to a two-dimensional square lattice with a 4-cell neighbor- 
hood, (b) The topology changes from a one-dimensional array to a two-dimensional 
square latticeto a 8-cell neighborhood, (c) Changing from a 4-cell neighborhood 
to a 8-cell neighborhood in a two-dimensional square lattice, (d) Graphical evi- 
dence that the Gini coefficient decreases in the transitions described in (a), (b), and 
(c): one-dimensional array (circle points, upper line); two-dimensional lattice with 
4-cell neighborhood (square points, middle line); 8-cell neighborhood (triangle down 
points, lower line). Fixed r = 10. 

In this paper we have characterized in detail a deterministic economic model 
implemented on a two-dimensional square lattice with periodic boundary con- 
ditions. The system depends only on two parameters; the parameter a brings 
information about the local environmental pressure while r expresses the self- 
growth capacity of each agent. By varying the values of these parameters, the 
model can exhibit asymptotic BG or Pareto distributions of wealth. The strik- 
ing fact is that only a change in either one of these parameters is enough to 
make the system undergo from one type of statistical behavior to a different 
one. We have called this process Transition Strategy I. 

Another situation that can produce a transition from Pareto to BG behavior in 
this system has been found in this work. We have called it Transition Strategy 
II. In this scenario, a variation of the values of the parameters is not necessary; 
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the transition in the type of asymptotic statistical behavior is a consequence 
of the rearrangement of the neighborhoods of each agent. In particular, it is 
observed that an increase in the number of neighbors of each agent drives 
the system toward asymptotic states where the wealth is distributed more 
equitatively, that is, with a lower value of the Gini coefficient. The specific 
regions of the parameter space where this transition can take place have been 
delimited. 

To our knowledge this is the first time in the field of economic systems where 
this kind of transitions is reported. We hope that this work can help to en- 
lighten the mechanisms that operate in the complex world of economical trans- 
actions among individuals, companies, countries or other economic entities. 
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